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On Selberg’s approximation to the twin prime problem 

R. Balasubramanian and Priyamvad Srivastav 


Abstract 

In his Classical approximation to the Twin prime problem, Selberg proved that for x sufficiently 
large, there is an n G (x, 2x) such that 2^0) _(_ < X with A = 14, where fi(n) is the number 

of prime factors of n counted with multiplicity. This enabled him to show that for infinitely many 
n, n{n + 2) has atmost 5 prime factors, with one having atmost 2 and the other having atmost 3 
prime factors. By adopting Selberg’s approach and using a refinement suggested by Selberg, we 
improve this value of A to about A = 12.59. 


1 Introduction 

The Twin prime conjecture is one of the oldest unsolved problems in Nnmber theory. A statement 
very simple to understand has eluded numerous attempts by the most adept of mathematicians. In 
fact, we seem nowhere close to settling this problem. The best known approximation to this problem 
is due to J.R. Chen [Che] in 1973, which states that there are infinitely many primes p, for which 
p + 2 has atmost two prime factors. 

In the last few years, there have been major developments on a related problem, called the Bounded 
gaps problem. This problem asks whether the quantity H^n = lim inf —Pn) is finite and provide 

n^oo 

an upper bound for the same. In 2006, Goldston, Pintz and Yildirim(GPY) in [GPYl] proved 

liminf ^"+^~^" =0 
n^oo log Pn 

thus settling a long standing conjecture. More interesting than their result was the method they 
used to arrive at this result. This method, known as the GPY sieve method, caught the attention 
of many experts in the held. In the same year 2006, Y. Motohashi and J. Pintz [GPY2] , intro¬ 
duced a smoothed version of the GPY sieve. Their method just fell short of proving the bounded 
gaps problem. It was hnally in 2013 that Zhang [Zha] showed Hi < 7 x lO”^, marking an important 
breakthrough in the subject. This bound was subsequently lowered to about 4680 by the Polymath 
8a [Pol8a] project. In his paper, Zhang used the GPY sieve, alongwith a modihed version of the 
Bombieri-Vinogradov Theorem. It was the modihed version of the Bombieri Vinogradov Theorem, 
which allowed him to deduce his results. 

Recently, in his preprint titled ‘small gaps between primes’, James Maynard [May] introduced 
a rehnement to the GPY method by using multidimensional sieve weights. Under this rehnement, 
Maynard shows that Hi < 600 and that Hm 'C This result was also proved independently by 

Tao. In his paper, Maynard has only used the Bombieri-Vinogradov Theorem and the paper does not 
incorporate any of the technology involved in the Zhang’s paper. The Polymath 8b [Pol8b] project, 
which extends Maynard’s methods, has successfully brought down Hi to about Hi < 246. Under the 
Generalised Elliott-Halbertsam (GEH) conjecture. Polymath 8b shows that Hi < 6. 



The use of multidimensional sieve weights actually dates back to Atle Selberg, who suggested their 
use in his approximations to the Twin prime problem [Sel, page 240]. In this unpublished manuscript, 
Selberg considers the sum 


E 

x<n<2x 
n=—l(mod 6) 



‘^(n) _j_ 2 ^{^^+ 2 ) 

A 


( V 


Ad 


( 1 . 1 ) 


where z = He succeeded in showing that (1.1) is positive for any A > 14. When this sum is 

positive, it means that there is an x < n < 2x such that 2^^”^ + < 14, which implies that 

n{n + 2) G P 5 . Since this happens for all x, there are infinitely many n for which n{n + 2) G P 5 . 

Note that if A can somehow be brought below 12 , one could show the existence of infinitely many 
n for which n(n + 2) G P 4 , since 12 = 2^ + 2^. 


By adjusting the sieve weights suitably, Gerd Hofmeister(unpublished) was able to bring down A 
to “about 13”, though Selberg is not clear about what exactly “about 13” means. 

Further, in order to improve the value of A, Selberg [Sel, page 245] suggested two-dimensional 
sieve weights of the form 


'y ^ ^di,d2 

di\n 

\d2\n+2 


( 1 . 2 ) 


In this paper, we adopt Selberg’s approach and improve the value of A by using the weights as 
given in (1.2). However, this approach falls somewhat short of bringing A below 12. 

We consider the sum 


E 

nr^X 

n=i;o(niod W) 



2^(71) _ j _ 2 ^(^+^) 

A 


/ 




Y AA,<i 2 


di\n 
\d2\n+h 




Si - S2/X 


(1.3) 


where h is an even number and uq is chosen such that {n{n + h),W) = 1. Following Maynard [May], 
we have defined 

D = log log log X and W = n P (1-4) 

p<D 


The sieve weights A^^^dj supported on 


S{z) = {{di,d 2 ) ■ fi^{did 2 W) = 1,max{d^^^(i 2 ,< z} 


(1.5) 


In this paper, we show that 
Theorem 1.1. (1.3) is positive for any A > 12.59. 


The paper consists of three major sections. In the first section, we introduce the Perron’s formula 
and use it to provide asymptotic estimates for partial sums of a certain class of arithmetic functions. 

In the second section, we deal with partial divisor sums like T{n) and in arithmetic progres¬ 
sions. This is because the quantity S 2 in (1.3) turns into a combination of sums of the function 
in arithmetic progressions. We write expressions for such divisor sums in arithmetic progressions in 
terms of suitable Kloostermann and Exponential sums. These results are essentially due to Selberg 










in his unpublished manuscript [Sel, page 241], The Weil’s bound for exponential sums allows us to 
obtain asymptotic formulae for the sum 


E 




n<x 

n=a(mod m) 


with a good error term. 

In the final section, we obtain asymptotic expressions for Si and S 2 (See 4.25 and 4.47) to compute 
a suitable value for A. In the computations, Corollary 2.14 from Section 1 is frequently invoked as 
it provides asymptotic formula for just the sums we encounter. Moreover, these calculations and 
estimates follow many ideas from Maynard’s preprint [May]. In order to optimise the choice of sieve 
weights, we rely on a sage program. (See 4.55). 


Notation 

Throughout this paper, h remains a fixed even number. The notation n ~ x means that n < x < 2x. 
The symbol p is reserved for a prime number. The numbers e, x and z will always be positive real 
numbers with z < x. Many times, we shall assume that x is sufficiently large and that e is sufficiently 
small. We write (a, b) and [a, b] for the GCD and LCM of positive integers a, b respectively. In many 
places, particularly Section 3, we write /(a,6)(or f[a,b]) to denote /((a,6))(or f{[a,b]) ) to simplify 
notation. For arithmetic functions / and g, f * g denotes the dirichlet convolution of / and g. The 
symbols 'O', ' denote the usual big oh notation and 'o' denotes the little oh notation. In many 

places, particularly Section 3, the 0-constants depend on e or the functions P, Qi or Q 2 , which 
will all be bounded functions. For the sake of simplicity, we do not specify the dependence of these 
0-constants. 

The following table has a list of arithmetic functions. All of these functions barring fl(n) are 
multiplicative. In some cases, we shall define tthe function on primes as we are only concerned with 
their value at squarefree integers. The functions /, /i, gi, h, hi and /i 2 appear only in Section 3. 


Function Description 


h 

id 

T 

n 

g 

f 

fi,gi 

h 

hi,h2 


The Moebius function 

Defined by id{n) = n 

Euler Totient function 

The Divisor function 

The additive function given by D(p*’) = k 

Defined by g{p) = 

Defined by f{p) = pl2 

Satisfy / = /i * 1 and g = gi*l 

h{p) = 1 - 3/(p -b 2) 

hi{p) = 1 — 3/p + 2/p‘^ and h 2 {p) = 1 — 2/p+ 2/p^ 
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2 Preliminary Results 

In this section, we review some of the well known results for partial sums of arithmetic functions. We 
define a class of arithmetic functions and give asymptotic estimates for their partial sums. We state 
a well known result 

Lemma 2.1. 

^!2ir = i„g:, + o(i) 

p 

p<x 

Proof. For a proof, see [Vau, Pg 50] □ 

This leads to an immediate Corollary. 

Corollary 2.2. Let P be a continuously diferentiahle function on [0, p] and let x^lf’ < z < x. Then 
we have 

log X 
log z 

Q 

Proof. We apply partial summation to the function . Since we know from Lemma 

2.1 that Y1 = logx + 0(1), we get 

p<x 


p 


Viog^y 


= /^ (Si) 


X 



f + o(i) + o(/V 



d 



Note that the main 
completes the proof. 


term above is 


clearly (logz) 



P(t) dt and 


the error term is 


0 ( 1 ). 


This 

□ 


We now state the following Lemma. 


Lemma 2.3. Let c > 0 and x > 0 be a real number which is not an integer. Then for any T > 0, we 
have 

I 


1 j-c+iT s 

^ / ^ds- 6 {x] 

ZTTZ Jc—i'j' S 


< X min < 1, 


T| logx| 


where 


5{x) = 


0 i/ 0 < X < 1 
1 if X > 1 


Proof. For a proof, see [Iwa, Proposition 5.54, Pg 151] 

This leads to an immediate Corollary, namely the Perron’s formula. 


□ 


Proposition 2.4 (Perron’s formula). Let F{s) he the Dirichlet series of f which converges absolutely 
for a > (Ja > 0. Let c > aaund let x > 0 be a non-integer. Then 




n<x 


1 

27ri 


pc-\-iT 

Ic-iT 


r(.)2-* + o(x'=f;hMmi„p 


n=l 


1 


’rllog^l 
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Proof. We have 


1 pc+iT s 1 fC+iT ( °° 

— / F{s)—ds = — / y 

2m J^_iT s 2m J^_iT 



Since the contour formed by the line joining c—iT and c+iT is compact, it follows that F{s) converges 
uniformly here. We can therefore swap the order of integration and summation to obtain 


hi 


c+iT s 1 

F{s)—ds = — 
c-iT s 2m 


n=l 


pc-\-iT 

1 


Jc-iT 

s 

\n) 


By Lemma 2.3, we have 


pc+iT 

2m s \n 


1 /X\s 


ds = 


l + 0(f^min{l,^l4^}) 


ds 


n < X 
n > X 


From the above two relations, the result follows. 


□ 


We define a class of multiplicative arithmetic functions and prove asymptotic formula for partial 
sums of the same. These are the type of functions we shall be encountering in the main computations 
of Section 3. 

Definition 2.5. For an integer A; > 1, we define to be the set of all multiplicative arithmetic 
functions / which are supported on the squarefree integers and the Dirichlet series F{s) of / is of the 
form F{s) = C^(l + s)G{s), where G{s) is given by an absolutely convergent series in cr > —<5, for 
some (5 > 0 . 


Proposition 2.6. Let / G be an arithmetic function with Dirichlet series F{s) = f^{s + l)G(s). 
Then for any non-integer x > 1, 

y fin) = + O x) 

n<x 

Note that here G{s) = 0 • 

Proof. From (13.10) of [Ivic, Pg 353], we have for any k > 1, 

y Tfc(n) = xPfc(logx)+ (2.1) 

n<x 


where 0 < 0 *, < 1 , = 1 * 1 * • • • * 1 and Pfc(log x) = i?eSs=iC^(s)^y is a polynomial of degree k — 1 


k times 
{\ogx)^~^ 


in logx, with leading term ^ • By partial summation, we obtain 


E 

n<x 


Tkin) 


n 


= 1 + 0 


log X 


(log a:) 
kl 


( 2 . 2 ) 


Since the Dirichlet series of Tk/id is Q^{s + 1), we have / = Tk/id * g. Therefore, by the convolution 
method, 

(log 2)1= 


E /(") = E «(“)^ = E i'f") E ^ = E <'(“) (1 + 0 

n<x ab<x a<x b<xla a<x 


1 

log x/a 


kl 


(logx)^ ^ ^ ^ loga^ ^ 


kl 


y 5(a) ( 1 


a<x 


logx 


(2.3) 


+ 0 


k-l 
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Since G{s) is well defined at s = 0, it follows that ^ g{a) = G{0) + A*{t), where A*{x) = 5(^)- 

a<t n>x 

Since G{s) = ^ converges absolutely for a > — 6 , it follows that A*{t) <C . Letting Q{x) = 


(1 — x)^, we have 


'^9{n)Q 


n<x 


log n 
logx 


/ 


= / Q 


logt 


logx 

= AQ{0) + Oix-^) + O 
= G{0)Q{0) + 0{x-^) 


diGiO)-A*it)) 


Q' 


, ( logt 


11 fi+'^logx^ \loga: 


dt 


□ 


Substituting the above expression into (2.3), we obtain the desired result. 

Remark 2.7. Actually, the asymptotic formula for /(^) ™ previous Proposition is of the 

n<x 

form CQ{logx) + 0{x~^), where C is the appropriate constant, Q is a polynomial and 0 < 0 < 1. 
Since the expression given in Proposition 2.6 suffices for our purposes, we avoid writing the main 
term as a polynomial and instead write with an error O 

We now redefine the constant G(0) occuring in the above Proposition. 

Definition 2.8. Let / G flfc and let / = /*!. Then for any positive integer m, we dehne 


p|m 


(2.4) 


Proposition 2.6 leads to the following Corollary, which gives us the asymptotic formula for partial 
sums of functions in 

Proposition 2.9. Let f ^ / = / * 1. Then for any positive integer m and {d, m) = 1, we 

have 

f{d) c(m,f) ( !og(iy 


E /( 


n] = 


n<z 

n=0{mod d) 

Proof. Firstly, we note that 


fid) 


m 


k\ 


1 - 


logz) 


+ 0(log^-i z) 




n<z 

n=0(mod d) 
(n,m)=l 


nKzjd 

(n,dm)=l 


We now apply Proposition 2.6 to the function f {n)l(^n,dm)=i with x = zfd. Note that the Dirichlet 
series of this function is C^(l + s)Gdmis), where Gdmis) = G{s) 0 therefore 

obtain, 

log'' z/d 


p\dm 


E = fid)^^^Gdm{0) + O (log"-i z/d) 


n<z 
n=0(mod d) 
(n,m)=l 


Writing 


o*,(0) = 0(0) n (1 + 0(1 - + /W) 

j[a) ra 


p\dra 




we obtain the desired result. 


□ 
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The next Theorem gives us the asymptotic formula for functions in accompanied by a smooth¬ 
ing function. 

Theorem 2.10. Let f G Qk and let f = f * 1. Suppose 0 < wi < W 2 ^ for some p > 0. Let 
d < W 2 and let P be a continuously differentiable function on [0,p]. Then for any positive integer m, 


log W2 
log i 


VogzJ \ m ) f(d)ik-l)\ J 

max|l^,l2£d) 

^ log z ’ log z J 


tt;i<n<R;2 
n=0{mod d) 
(n,m)=l 


^ogd \ 

log zj 


k-l 


dt 


+ 0(log^-^ z) 

Proof. Let fm,d be dehned as in the previous Proposition 2.9. The required sum to be estimated then 

logn' 


IS 




n<z 


log z 


We have seen in Proposition 2.9 that for any x > 0 

X] fmA^) = M{x) + E{x) 


n<x 


where M{x) is the main term and E(x) is the error term. We also know that E{x) <C log^ ^ x and 
that M{x) is a differentiable function on We apply partial summation to get 


^ fm,dA)P 

wi<n<W 2 


log n 
log z 


Ul2 


= / P 


logt 

log z 


diM{t)+E{t))dt 


Wi 


The main term above is 


and the error term is 


W 2 


M'{t)P 


log^ 

logz 


dt 


Wi 


W 2 


O \E{w2)\ + \Eiwi)\ + 


' \m\ 

tlogz 


dt 


Wl 


Since E{t) <C log^ ^ t, it follows that the error term above is O (log^ ^ z). Substituting the value of 
the main term M{t) from Proposition 2.9, the main term above is 


1 fjd) f y?(m) log z 

k\J{d) V nfi 


u '^2 

\ ..X f flogt logd\ 

) J Vlog^ log2:/ 


fe-i 


P 


maxjuii ,d} 


logt\ k dt 


log Z J t log Z 


The change of variable t ^ yields us the desired main term. 

log z 

Definition 2.11. For any real number s, we define 


One easily sees that 


, , . J 1 2s 3(1 - s) 

?](s) = mm <; 1 - - - - 


1 2s f 

fcii .f»>5 

2 “5 


□ 


(2.5) 


r?(s) = < 
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We state the next Lemma without proof. 


Lemma 2.12. Let di, d 2 be positive integers. Let z < x, S{z) he as given in (1.5) and W he as given 
in (1.4)- Then 

{di,d 2 ) G 5 ( 2 ;) p?{did 2 W) = 1 and vii&yi{did^^d 2 } < 2 ; 

/ log (j; 

pL^{did 2 W) = 1 and di < 2 :^V 
Definition 2.13. Let si,S 2 G [0,1]. We define the region 


Tsus2 = \ {xi,X2) G : Si < Xi, xi + ^ < 1, ^ + X 2 < 1 


( 2 . 6 ) 


We shall denote the region Tq^q by T. 


The next Corollary is the key result of this section. In our computations later on, we shall be 
invoking this lemma quite frequently. 

Corollary 2.14. Let fi G Llki and f 2 G 11^2 fox positive integers ki and k 2 . Let f = fi * f 2 and 
let P : T ^ M. be a function differentiable in each variable. Let di,d 2 be positive integers and let 
~ ^ ^ ~ triw^ogz ^ Then with S{z) as defined in (1.5), we have 


E h(h)Hh)P (kik. kik) ^ Sff) AW y) 

^ ylog^ log^y f{di) /(da) 


li=f){mod di) 
(hM^S{z) 




V Tji,S 2 

+ ©(log'^l+^^-l 


(fci - 1)! (fca - 1)! 




Proof. First of all, let us rephrase the conditions h = 0(di), I 2 = 0 (^ 2 ) and {li,l 2 ) G S{z). We know 
from Lemma 2.12 that 

{li,l 2 ) G S{z) ^ p^ihhW) = 1, and h < zWiog.) 

From now on, we let L = ^ \ for i = 1,2. We can now write the given summation as 

log 2 ; 

E = E E = E E 

q=0(di) h=0{di) l2=0(d2) q=0(di) l2=0(d2) 

l2=0{d2) (h,d 2 )eS{z) {l 2 ,h)£S(z) l 2 <z^Oi) 

{h,l 2 )&S{z) (li,d 2 W) = l {l 2 ,hW) = l 

This given summation above is set in such a way that one can directly apply Theorem 2.10 to the 
inner sum. We therefore have 

E E Aft) E 


li=0{di) 

(ll,l 2 )eS{z) 


/2(d2) 


Zi= 0 (iii) l 2 = 0 (d 2 ) 

Zl<2’>ft2) l2<Z^Ol) 

{h,d2W) = l {l2,hW) = l 


log z ’ log Z J 


^ 2 (^ 2 ) 

Zi< 2 ' 7 (s 2 ) 

l,h42W) = l 


E Aft)(®^)‘k(.ri../,) 


log ^1 

log 2 


P 


(2.7) 


log 2 :’ J (^2-1)! 


52 


+ 0 I (logz)'=^-i J]/l(/l) 

ll<Z 
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We can now write 


iwh,f2)= n (i-^) 'Mp) 


p\Wh 


Pih)J pj 


1 \ 


k 2 


n 1-- mp ) 


AWJ 2 ) 


Mh)J f2ih) 

Therefore, substituting this expression for c(W/i,/ 2 ), the main term in (2.7) may be written as 


/2(rf2) / V?(PE)log2: 

/2(d2) V ^ 


k 2 


iw,f 2 ) 


h=0{di) 

{li,d2W)=l 


hih) 

f2{ll) 




p 


Vlogz’ (^ 2 - 1 )! ^ 


S2 




if 2 /f 2 ){d 2 ) ((p{W)logZ 


(/2//2)(d2) V W 


ki+k 2 


c(W,f 2 )ciW,^ 


V(s 2 )v(tl) 


I I P{tl,t 2 ) 


{h - si)P ^ A - S2A ^ 
{h - 1)! {k 2 - 1)! 


dt 2 dti 


Since 


(/i//i)(di) 


(/2//2)(d2) flA)f 2 id 2 ) 
{f 2 /l 2 ){d 2 ) f(di) f{d 2 ) 


and c(W, f 2 )c {W,=J= c{W, h * f 2 ) = c{W, f) 


we obtain the desired main term. Moreover, since fi G it follows that the error term in (2.7) is 
0 (log*^i+^ 2 -i completes the proof. □ 


3 Estimates on Divisor sums 

The computation of S 2 transforms into certain divisor sums in arithmetic progressions. So, in this 
section we shall obtain asymptotic estimates for the sum 

Dm,aA = ( 3 - 1 ) 

n<x 

n=a(mod m) 

where m is an even squarefree integer and {a,m) = 1. To do so, we give an expression for Dm,a{^) 
in terms of Exponential and Kloosterman sums (See Proposition 3.8) and later use the Weil’s bound. 
Using the convolution method, we can the use this expression of Dm,a{^) asymptotic formula 

for the sum 

c^{n) 

n<x 

n=a(mod m) 

These results can be found in Selberg’s manuscript [Sel, Pg 234-237] 

Definition 3.1. Eor any positive integer m, define 

^m{x) = ^ r(n) 

n<x 
(n,m)=l 



(3.2) 
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The next Proposition gives us an expression for Am{x). 

Proposition 3.2. The following estimate holds for any squarefree positive integer m < x 


Am{x) = X- 




nogx + c + 2 ^^^ j +0, 

\ p|"i / 


Proof. Define a function by 


Tm{n) = 


T{n) if (n, m) = 1 
0 otherwise 


Then the Dirichlet series for r™ is 


Fm{s) — 


r(n) 


(ri.,m)=l 




pfm 


- n ^ ~c^{s)Gra{s) ( 3 . 3 ) 


p\m 


Let gm be the function defined by the Dirichlet series Gm{s). Since C^('S) is the Dirichlet series of 
T(m), it follows from (3.3) that 

Xm = T =f= Qm 

Moreover, gm is given by 

9min) = T{a)h{b) 


a^b\m 

ab=n 


Using the convolution method, we therefore have 


E 


Xm[n) = YJ2'^\d) 
n<x n<x d\n 


d<x 

— — a 




d<x 


x(log X — log d + 27 — 1) 


d 


+ 0 


T/2 


= X 


(log x + 27 - 1 ) Y 


. d<X 


gmjd) \ 

, j 


— X 


dV 2 ^ 
gm{d) logd 


d<x 


d<x 


Now, 


We therefore have, 






Also, 


Therefore, 


E 


gmjd) logd 

d 


-G71) = - 2 G.„( 1 )E^ = - 2 ^E‘“" 

z_—/ p — rr) ‘^ 


p\m 




p\m 


E 

d<x 


gm (d) log d (d) log d 


d 


E 


d 


+0 E 


\d>x 


gmjd) logd ^ 


d 


^ ^ logP f x^(m) log X 




p\m 


p-1 


(3.4) 


ymyu) mg u ^ ^ j ^1/2 gmjd) j 




v« 2 iM = E( 7 )+o(y^) = EW+o(LM) (3.5) 

^ d w? \ ^ d / \ X J 

d<x \d>x / ^ ^ 


(3.6) 
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and the error term in (3.4) is 

«E E ^ E E (3-7) 

d<x d<x a,h\m a,b\m 

Substituting the relations (3.5), (3.6) and (3.7) back into (3.4), we obtain the desired result. □ 

Next, we give a relation of Divisor sums Dm,a{x) and Am{x) (See 3.1 and 3.2) in terms of Expo¬ 
nential and Kloostermann sums. 

Notation 3.3. By e(a), we shall mean exp (27riQ;) 


Definition 3.4 (Kloosterman sums). We define 

S{a,b,m)= ^ 

/ih=l(mod m) 

A sum of this type is called a Kloosterman sum. 


ah + bh 


m 


(3.8) 


We now state a famous result due to A. Weil. 
Theorem 3.5 (A. Weil). 


\S{a,b,m)\ < m^/^r(m)(a, 6, m)^/^ 


Definition 3.6. Eor |a| < -, define 


St{a) = e{na) 


l<n<t 


Lemma 3.7. 


|5i(«)l < 


2 a 


Proposition 3.8. Let m be any positive integer and {a,m) = 1. Then we have 


f r. , ^dt 


E 

n<x 

n=a(mod m) 




n<x 
(n,m)=l 



1 


(3.10) 
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Proof. First we rewrite the RHS of (3.9) in the form 

X r 

X] S{ac,d,m)St (—) S^/t (— 

The sum inside the integral above is 

d 

m 


dt 

t 


(3.11) 


S{ac,d,m)Sj 

' ach + dh 




s<x/t 

Recalling the fact that 


equation (3.12) becomes 


E E' 

-f<c,d<f hez^ 


E E E 

r<t y_2n<c<^ 


m 


E' 

r<t 

s<x/t 


rc + ds 


m 


(3.12) 


c{ah + t] 


m 


E 


-f<d<f 


d{h “h 5^ 


m 


E = 

6(mod m) 


m if X = 0(mod m) 
0 otherwise 


m 


Y. Y l = m2 ^ 1 

r<t r<t 

s<xlt ah=—r{raod m) s<x/t 

rs=a(mod m) ;i=-s(mod m) rs=a{mod m) 

The cardinality of {h G : ah = —r(mod m), h = —s(mod m)} is 1 since ha = —r(mod m) implies 
that h = —ra(mod m) and this forces a unique choice for /i(mod m). 

Therefore, (3.11) becomes 

/ \ 


E 1 


r<t 

s<x/t 


dt 


Y 


rs<x 

rs=a.(mod m) 


^ T(n) log 


\rs=a(mod m) j 

This proves (3.9). 

To prove the next part, we rewrite the RHS of (3.10) as 


n<x 

n=a(mod m) 




1 l-f<c,d<f 


^ S’(ac,0,m)5((i,0,m)S’i ^ (—) 

rn ^ \m J 


dt 

t 


(3.13) 


Again, the sum inside the integral above is 

S{ac,0,m)S{d,0,m)St(^-Y 


-f<C,d<f 


E E ^ 

— ^<c,d<^ h\,h2&Z^ 

E E f E ' 

S<xlt 


achi + dh 2 


m 


E' 

r<t 

\s<x/t 


rc + sd 


m 


/ 


c{ahi + r] 


m 


E 

m ^ m 

-^<d<^ 


d(/l 2 + s] 


m 


(3.14) 
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Again, recalling the fact that 


equation (3.14) becomes 



m if X = 0(mod m) 
0 otherwise 


r<t h\,h2&Z^ 
s<x/t ahi=—r{raod m) 
/i2=—s(mod m) 


= 1 

r<t 

s<x/t 

(rs^m)=l 


Note that |{(/ii,/i 2 ) G = —r(mod m), /12 = —s(mod m)}| is 1 since ahi = —r(mod m) 

implies that hi = —ra(mod m) for which there is only one solution. Similarly, there is a unique choice 
for /i 2 . Hence there is exactly one choice each for hi and /i 2 modulo m. 

Therefore, (3.13) becomes 


X 


k 


E 1 

r<t 

s<x/t 

(rs,m)=l 


dt 

t 


E 

rs<x 

(‘rs,m)=l 


n<x 

(n,m)=l 


This completes the proof of the Proposition 3.8. 


□ 


Lemma 3.9. Let m be a squarefree positive integer. Then 

S{a, 0, m)(5,0, m) 


S{a, b, m) — 


(p{m) 


< 2mVV( m){a, b, m)^/^ 


(3.15) 


When either of a or b is = 0{mod m), the above difference is 0. 


Proof. First, note that when either of a or 6 is = 0(mod m), say b = 0(mod m), then the quantity in 
the LHS of (3.15) is 

Qf n 'i 5(0,0,771)5(0,0,m) 
b{a, 0, m) - 


ip{m) 


= 0 


since 5(0,0, tti) = ip{m). 

When neither of a or 6 is divisible by m, we apply the trivial estimate. By Theorem 3.5 

|5(a, 6, r?i)| < ?77^/^T(77i)(a, 6, ?7i)^/^ 

Secondly, since m is squarefree, we have 

S{a,0,m) = dp(^—'^=iJ.{m) dfj.{d) = 


So, 

Therefore 


d 

d|(a,m) d\{a,m) 

|5(a,0, ?7i)| = ip{a,m) 


(3.16) 


5(a, 0, m)S{b, 0, m) 


qp(m) 

This completes the proof. 


(p{a,m)(f{b,m) ^ (p{a,b,m)(f[{a,m), {b,m)] ^ ^ 


ip{m) 

< (a, b, 777 )^/^ 771 ^/^ 


ip(m) 


□ 
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We shall make use of Propositions 3.2 and 3.8 in order to obtain an expression for Djn,a{x) with 
a decent error term. 


Theorem 3.10. Let m be a squarefree positive integer and let {a,m) = 1. Then for any e > 0, we 
have 


Dm,a{x) = I log X + C + 2 

p\m 


whenever m < ^. 

Proof. Consider the difference 


logp \ „ / x^/ 2 +€ 

p-l}^ ^ \ mV4 


Tim,aid) / 

(p[m) 


dt 

t 


(3.17) 


(3.18) 


From Proposition 3.8, the above expression (3.18) is 




^ / S’(ac,0,m)5(d,0,m)\ /c\ fd 

^ ( S{ac,d,m) --) St [-)S,/t ( - 


rf<c,d<f 


m 


m 


dt 

T 


(3.19) 


By Lemma 3.9, the quantity inside the integral in (3.19) is zero when either c = 0(mod m) or 
d = 0(mod m). So we may assume that neither of c, d is 0. Then we have 


S{ac, d, m) — 


5(ac, 0, m)S{d, 0, m) 


(p{m) 

Moreover, by Lemma 3.7, we have 


< 2m^/^T(m)(ac,d,m)^/^ < 2m^/^r( m)(c, d, 


St 


m 


m 

< —TT and 
“ 2 c 


m 


< 


Therefore, (3.19) is 


1 

< — 


^ 2 m^/^r(m)(c, d, 


m 


c,d^0 


4|c||d| 


m 

W\ 


dt 

t 


<2m'/V(m.)logi ("'4“)*'"" <2m‘/V(m)logg 

C-Ci OCX 


0<c,d<^ 


0<c,d<^ 


(3.20) 


Observe that 


E 

0<c,d<- 


cd 


/ ^ ^ ^ (d d^ 

g<f 0<c4<^ 9<f 0<c',d'< 

{c,d)=g 




E 0-^'^- 


E fj7 <':(3/2)log- 


m 




<C log^ X 

Therefore, from (3.20) and (3.21), we get 


(3.21) 


T{m) J t 


<C m^/^r(m) log^ x <C 
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We now let 


A{x) = Dm,a{x) and B{x) = 


ip{m 

l /2 + e 
m- 


'-^m (^) 


Let C{x) = A{x) — B{x). We want to show that C{x) <C ^^ 1/4 • Let 

5 = (3.22) 

It is then clear that 5 > m whenever m < Since dt <C it follows that 


px^5 

J x—S 


-S t 


(3.23) 


Note that we have 




n<x 

(n,m)=l 


^=a(mod m) 


ip(m] 


Therefore, for any interval of length m, the jump in C{x) is utmost Now suppose that C{x) ^ 

^^ 1^4 for some x. Consider the neighbourhood Ix = {x — d,x + 6 ) of x. We then have for any t ^ Ix 
that 


C{t) > 


d/ 2 +e 


m 


1/4 


+ O —x 


m 


Under this choice of d in (3.22), one immediately sees that C(t) !:$> ■ It then follows that 


L 


x-^S 

-s t 


C(t) X" 

^ ^ dt:^ 


l/2+e i-x+S , f X + 5 


m 


1/4 


(M 

J x—S t 


m 


1/4 


log 


X — 6 


> 5- 


,- 1 / 2 +e 


m 


1/4 


= 


This is a contradiction to (3.23) and the proof is complete. 


□ 


Theorem 3.11. Let m = lUm^ he an even squarefree integer, where m' is odd and {W,m') = 1. 
Suppose (a, m) = 1. Then for any e > 0 


V = , v(tr) C{W) 

^ W 2 g{m') 


( 


n<x 

n=a{mod m) 


log X+c+2 


V 


p\Wm.‘ 

p>2 


logp 


+ 0 . 


r.l/2+e 


m 


1/4 


(3.24) 


n (p ^ ^ P ^ ) 

—7 -- and g(p) = -. 

pip — 2) p — 2 

p\W ^ 

Proof. We shall make use of the Convolution method. We write 

n 
d 


d\n 


Then the Dirichlet series of an is given by 

(1 - p -^)2 


ki(s)= n 


(1 - 2p-U 

p \ J. / p 




1 


p^{p^ - 2 ) 


= Ci2s)Bis) 


It then follows that > 0, for all n. 

Note here that B(s) is convergent for a > 1/3. Observe that A(s) has a simple pole at s = 1 / 2 . This 
means that 


E 

n<x 


ni/2 


O(logx) 


(3.25) 
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We have 


n<x 

n=a(mod m) 


T 


n 


E E E«-(5)-E- E -u 

n<x d\n d<x n<x 

n=a(mod m) (d,m)=l n=a(mod m) 

n=0(mod d) 

^ X] ^ = X] 

d<x ^ 

(c/,m)=l n=a^(mod m) (d,'m)=l 




E 

d<x 

(d,m)=l 


ad (fj-m) 
d m? 


[^log X - log d + c + 2 ^ 

V 


+ o. 


2;l/2+e 


E 


d<x 


From (3.25), the error term above becomes O logx^ = O The main term is 


(p{m) 




\ p\m J d<x 

(d,m)=l 


ad ‘fim) ad, , 

-X -^ 2^ — logd 




d<x 

{d,m)=l 


(3.26) 


Now, 


Moreover, since 


E ad _ I (Q ' ^ 

~d ~ ^ ~d^ 

d<x (d,m)=l 

[d,m)=l 


TT (p-l)% o/!ggf 


X I ^^p{p-2) 
p\m 


1 


—=TT (1+■ 

d^ ill pHp^-2) 

{d,m)=l p\m 


we can take logarithmic derivatives to obtain 


(3.27) 


d" 

{d,m)=l 

Snbstitnte s = 1 to obtain 


E gi°g<i=-n(i+ p.(p(_2) ) (~^E 


p\m 


logp 


p\m 


(P‘ - 1)(P* - 2) 


w ^iogd = TTk:ilL 

E^ d p{p — 2) 


(d,m)=l 


p|m ' 


2 ‘ ^ logp \ _-[-[- (p - 1 )^ 

^(p-l)(p-2) I iiMp-2) 

pfm / p\m 




C 0 - 2 E 


logp 


V 


p|m 

p>2 


(p-l)(p- 2 ) 


Snbstitnting the eqnations (3.27) and (3.28) into (3.26), we obtain 

( 


y = 

E^ -I--I- p{p — 2) 

n<x p\m 

n=a(mod m) 


log X + c + 2 'y 


\ 


logp 

I ^ 
p|m 

p>2 / 


+ 0 . 


„l/2+e 


m 


1/4 


/ 

(3.28) 


The proof is completed with the following observation 


^ _ TT (P- 1) TT (P- 1)^ _ P(^) 1 TT (P - 1)^ _ P(^) C(1T) 

li p2 lip(p_2) VF^ gim') -l■-l■p(p —2) VF^ gim') 

p\m p\m p\W 




□ 
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We have an immediate Corollary from the above Theorem. 


Corollary 3.12. Let m = Wm' he an even squarefree integer, where m' is odd and {W,m') 
Suppose {a,m) = 1. Then for any e > 0 and m < , we have 


= 1 . 




Lp{W) c{W) 
g(m') 


( 


n=a(mod m) 

where n ~ x means that x < n <2x. 


log X + c' + 2 


V 


p\Wm' 

p >2 


logp 


\ 


+ Oe 


^ 1/2 + 6 


m 


1/4 


4 The Modified Sieve 


We now use two-dimensional sieve weights to Selberg’s approximation of the twin prime problem and 
reduce the value of A. We shall first provide asymptotic expressions for Si and S 2 in Proposition 
4.3. The later part of this section will be devoted to proving this Proposition. Once these proofs are 
completed, we make use of a sage program to optimise the sieve weights and determine a suitable 
value of A. Following Maynard [May], we set 

D = log log log X and W = n p 

p<D 

This is done in order to avoid complications in our calculations. 


Let h be an even number and vq be chosen such that (n, IT) = (n + h, W) = 1. We consider the 


sum 


where 


and 


Let 


E L 


2^(71) _ j _ 2 ^(^+^) 

A 


/ 


n=i;o(mod W) 


\ 




d2 




di\n 
\d2\n+h 


Y 

di\n 


= Si- S2/X 


d2 


n=vo(mod tV) / 



f \ 

S2 = 

Y 

Tl'^X 

n=i;o(niod W) 

di\n 

\d2U-l-fc / 


(4.1) 


(4.2) 


(4.3) 


(4.4) 


— {idi,d 2 ) ■ fi‘^{did 2 W) = 1, maxjdid/^,< z} 

This is going to be the support of the sieve weights \di,d 2 - reason for this choice will soon be 
clear. We will choose 

2 = ^ 45 ^ 


Notation 4.1. Until the evaluations of and S 2 are complete, the set {i,j} will be a permutation 
of the set {1,2}. So, if we write Sij \ {di,lj), it means that both the conditions S 12 | {di,l 2 ) and 
•521 I {d2,h) hold. Moreover, if we write di = 0(mod k), it will mean that both the relations 
di = 0(mod li) and d 2 = 0(mod I 2 ) hold. 
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Instead of making a choice for \d^_,d 2 directly, we define the quantity 


D 


ri,r2 


f(ri)f(r2) ^ 

di=0(mod 


^di ,d2 

, fidi)f{d2) 


and make a choice for Zlri,r2- Note that -Dri,r2 vanishes whenever (ri,r2) is not in S{z). 
Here, / is defined multiplicatively by 


(4.6) 


f{p) = p/2 


(4.7) 


We set 


^ri,r2 — ^ 


piri)iJ,{r2)P 


0 


log ri log r 2 
logz ’ logz 


(ri,r 2 ) G S{z) 
otherwise 


(4.8) 


where P : T —)■ M is a bounded and symmetric differentiable function, with T as defined in (2.6). 

The reason for choosing P to be symmetric is that we want the Xdi,d 2 '’^ be symmetric and this 
is equivalent to the being symmetric. Morevover, choosing Pr-i,r2 is equivalent to choosing 

Xdi,d 2 due to the relation 


^di,d2 

f{di)f{d2) 


p{di)p(d2) ^ 

/i=0(mod di) 


p{h)p{l2) p 

f{h)f{h) 


p{di)n{d2) Y1 


/i=0(mod di) 

ihM^siz) 




I' logh log^2 

\ log z ’ log z 


(4.9) 


Notation 4.2. We set 


B = 


^p{W) log z 
W 


(4.10) 


The following Proposition gives us the asymptotic expressions for S'! and S 2 in terms of the 
function P. 

Proposition 4.3. Let Xd^^d 2 be as in (4-9) and let rj he as defined in Definition 2.5. Then with 

z = 

S, = ^{l + o{l))B^Ri{P) 

9 t 

S2 = ^{1 + o{1))B‘R,_{P) 

where 


Here 


RiiP) = jj Q\{spS2)ds2 dsi 

T 

R 2 iP) = jj (si(3 - si)Q 2 (s 1)'S2)+4siQi(si,S2)Q2(si,S2))ds2 dsi 

T 


(4.11) 


v{si) 


Qi{si,S2) = j j P{ti,t 2 )dt 2 dti and Q 2 isi,S 2 )= J P{si,t 2 ) 


dt2 


52 


(4.12) 
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The remainder of this section will be devoted to proving Proposition 4.3. Before we begin, we 
define two new quantities which will arise when computing Si and 82 - Let 

^di,d 2 


Brur 2 = ‘p{ri)^{r 2 ) 


di=0(mod Vi) 


did 2 


Crur 2 = fi{ri)gi{r 2 ) ^ 


^di,d 2 


f{di)g{d 2 ) 


di=0{mod Ti) 

where / = /i * 1 is as defined in (4.7) and g, gi are defined by 

^ and gi*l= g 


(4.13) 


(4.14) 


Lemma 4.4. The following relations hold for any (ri,r2) G S{z) 

a) |Ari,r2l log^ -2; 

b) Br^ r 2 = Ki{W)B‘^gL{ri)gL{r 2 )h{ri)h{r 2 ) Qi () +0(logz) 

\ log z log z / 

where Qi is defined as in (4-12). Here 

ki(1 T) = ]^ Tl --'j fl + -'j and h{p) = 1 - 3/{p + 2 ) 

p\W ^ PJ \ PJ 

c) C'ri.ra = B /i(ri)/x(r 2 )hi(ri)/i 2 (r 2 ) Q 2 +0(1) 

where Q 2 is as defined in (4-12). Here 

/ii(p) = 1 — 3/p + 2/p^ and /i2(p) = 1 - 2/p + 2/p^ 


Proof. 

a) From (4.9), we have 


A 


ri,r 2 


T{ri)p{r2)f{ri)f{r2) ^ 


li=0{ri) 

{h,h)&S{z) 


f{h)f{l2) 


f ^ogh log/2\ 

V log z ’ log z / 


Applying Corollary 2.14 to the RHS of the above expression and writing s, = we obtain 
Ari.ra = c{W, f * f)B^ ^2)(tl “ Sl)(t2 “ S 2 ) ^^2 + 0(log^ ^) 

p\r\r2 Ti,s2 

< log^ Z 

Here c{W, f * f) is as defined in (2.4) and tends to 1 as x ^ 00. 


b) From (4.13) and (4.9), we have 


Bri,r 2 

(p{ri)(pir2) 


E 

di=0(mod Ti) 


^di,d 2 

did 2 


E 

di=0(mod 


n) 


p{di)p{d2) 

T{dl)T{d2) 


E 

/j=0(mod di) 

{hM^s{z) 


ilMlMp 

fih)f{l2) 


E 

/i=0(mod Ti) 
ihMeSiz) 


f{h)f{l2) 


f log h log I 2 \ p{di)p{d 2 ) 
Vlogz’logz/ ^ r(di)r(d2) 


Zi=0(mod Ti) ^ ^ 

{h,h)&S{z) 


f ^ogh log I 2 \ 

V log 2: ’ log 2 / 


/logji log I 2 
\ log 2 ’ log 2 
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Applying Corollary 2.14 to the last equality above functions iJ?/id, iJ?/id G , we obtain the 
desired result. 


c) Again from (4.13) and (4.9), we have 


a 


ri,r 2 


E 


A, 


di,d 2 


fi{ri)giir2) .f{di)gid2) 

ai=0(mod Vi) 


, , . Id{d2)f{d2) iJ?{h)g?{l2) { logh log/2 \ 

d.Mdr, U^dd, Vlog.’log.; 

(h,l 2 )&S{z) 

g^ih)fj-‘^{l2) flogh log/2^ ,^^ld{d2)f{d2) 

E — ji :.—E >*(*)■ 


hk 

fiin) n{r2)ip{r2) 


/i=0(mod Vi) 
{hMeSiz) 


E 


i^Hh) 


ri\di\li 


g{d 2 ) 


fin) gir2) ^{h) 

/ 2 = 0 (mod r 2 ) 
iri,l2)SS(z) 


In the last equality in (4.15), we note from Lemma 2.12 that 

( log ri \ 

_ log^ 7 

By applying Theorem 2.10 to the sum in the last equality of (4.15) with the function G Hi, 
we shall obtain the desired result. 


□ 


We now begin with the evaluation of Si. 


From (4.2), we have 


4.1 Evaluation of Si 


Si 


A' 


E 

U'^X 

n=i;o(niod W) 


^di,d 2 

di\n 

\d 2 \n+h 


E 


E 


,(^2 I 2 


n~3; 

n=j;o(mod W) 


di li \n 
d 2 ,l 2 \'n+h 


^ ^ ^di,d2^h,l2 ^ ^ 


1 


{di,lj)=l 


n~x 

n=0(mod [dili]) 
n=—h{mod [d 2 ,L]) 
n=Do(niod W) 


(4.16) 


In the above sum, we have the conditions di,li \ n and ^ 2^2 | n + h. We first choose x large enough 
so that D > h. This ensures that rad{h) \ W. It then follows that {di,dj) = {di,lj) = 1 for i 7 ^ j 
and that {did 2 lil 2 ,W) = 1. This is because if there is a prime p dividing {di,dj) ( or {di,lj)) for 
i 7 ^ j, then p must divide both n and n + h and therefore p \ h. But since rad{h) \ W, it follows 
that p I W. This is a contradiction because the numbers di and Ij are all coprime to IT. Moreover, 
the conditions {di,d 2 ) = 1 and (/ 1 S 2 ) = 1 can be dropped because they are already included in the 
definition of S{z) (See (4.4)), the support of the sieve weights Xdi,d 2 - So, we are only left with the 
conditions {di,l 2 ) = {d 2 ,li) = 1- Since the numbers [di,/i], [ 42 )^ 2 ] and IT are pairwise coprime, the 
inner sum in the last equality of (4.16) becomes 

X 


IT[dl,/l] [42,^2] 


+ 0 ( 1 ) 
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Therefore, 


'S'l ^ ^ ^di,d2^h,l2 

{di,lj)=l 
= Ml + Ell 


W[di,li][d 2 , h] ^ W ldiMd 2 ’X] ^ ^ (155 

(4.17) 


The error term in (4.17) is 


{ X] I <C |S’(2;)plog®2 < 2;®/^log®z 


(4.18) 


. dl,d 2 


,(di,d 2 )eS(z) 


where we have used the bound |Arf^^rf 2 l ^ log^z from Lemma 4.4 and the fact that | 5 ( 2 :)| <C 


In the main term Mi, we can write 


1 

\di , /j] 


1 

ddi 


E 

^i I (^i i^i) 


v{n) 


This gives 




E 

{di,lj)=l 


Adj ,^2 ,L 

did2 I 1 I 2 


(p{ri)(p{r 2 ) 

ri\{di,li) 


To get rid of the conditions (di, ^ 2 ) = (<^ 2 , ^i) = 1, we multiply a factor of 


Y /^(Sl2)d(s2l) 
I (^* lb ) 


This idea was used by Maynard in his preprint [May]. 

Now, ri I (diji) and Sij \ {di,lj), for i / j. This implies that {ri,Sij) = {ri,Sji) = 1 because if 
there was prime p dividing and Sij, then p \ This is contradiction since [h^h) = 1- Also 

note that both and Sij divide di. Since {ri,Sij) = 1, it follows that riSij \ di. Similarly, we get 
TjSji I dj. Therefore it follows that both (risi 2 , r 2 S 2 i) and (riS 2 i, r 2 Si 2 ) are in S{z). Summarising 
this, we have 


n di,li= 0 {ri) 


^di,d 2 ^li,h 

did 2 I1I2 


Y /^(■S 12 )m('S 2 i) 

^ij I {di ,lj ) 





E 

k= 0 (riSji) 



(4.19) 


Remark 4.5. We need not write the conditions (risi 2 , 5 " 2 S 2 i) G S{z) and (riS 2 i, r 2 Si 2 ) G S{z) at 
every step because one of the two bracketed quantities in the last equality of (4.19) will vanish when 
any of these conditions do not hold (This is because the support of Xcp^d 2 is S{z)). 

Recalling the Definition of Rri,r 2 from (4.13), note that the two bracketed quantities in the last 
equality of (4.19) can be replaced by appropriate Rri,r 2 ’s i-e we obtain 


Ml 


X ^ /i^(ri) p^{r2) p{si2) p{s2i) 

^ 4 ^ ‘P(ri) <f(r2) ^ <f^(si2) P^(s2i) 


(4.20) 
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We would now like to get rid of the numbers S 12 and S 21 from the above summation. Since 
('512,521) G S{z), we have (.si 2 ,W) = {s 2 i-,W) = 1. Therefore, either Sij = 1 or Sij > D. The 
contribution to the main term Mi of the above sum (4.20) when atleast one of Sij is greater than D 
is 


<C 


< 


E 


1 


E 




^ I ^ ^ ^{n) ip{r2) 

(si 2 ,S 2 i)&S{z) {ri,r 2 )€S{z) 

S12>D 


X log^ Z 

W 




E 

ri,r 2 <z 


1 


<p{r)(p{r 2 ) 


< 


X log° 2 ; 

WD 


where we have used the estimate i?ri,r 2 ^ log^ -2^ from Lemma 4.4. 

We may therefore assume that S 12 = 521 = 1 with a cost of error term of 


Eu = 0 


X log® 2 ; 
WD 


Therefore, one can now write 


^ X fi^{r2) ^2 , n, , ^ 

- iv E ^(10 

Substituting the value of i?r-i,r 2 from Lemma 4.4, we obtain 


X 

W 


/i2(ri) /i2(r-2) / 2 .wid4,2 


^1 = 777 E E:VEEiM ^i(^)^^^^(n)/^^(r2)Q^ 


{ri,r 2 )&S{z) 


‘Pin) pir2) 


2 [ log ri log r 2 


log 2 ; ’ log 2 ; 


(4.21) 


(4.22) 


+ 0(log3 z) (4.23) 


Note that the inner sum in (4.23) is precisely the type of sum estimated in Corollary 2.14. Applying 
Corollary 2.14 with the functions h?/p, ^ jp G fli ( See Definition 2.5), we obtain 


Ml = — B® 

IT 


Qi( 5 i, S2) ds 2 dsi + O 


X log® 2 ; 
IT 


(4.24) 


with Qi dehned as in (4.12). 

When we choose 2 : = we have 


En «and 


This gives the following final expression for 5i 


where 


5i = -(l + o(l))i?® i?l(P) 


RiiP) = // Qlisi,s2)ds2dsi 


(4.25) 


( 4 . 26 ) 
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4.2 Evaluation of S 2 

We will use the same ideas here as we did while computing 5i. The conditions and observations given 
in the paragraph right after equation (4.16) will be implemented. 

From (4.3), we have 


52 


E (2^W + 2^(-+'*)) ^ 


ri'^x 

n=t;o(niod W) 


di ,/i \n 
d2,l2\n-\-h 


E E ( 2 ^W+ 2 ^(-+^)) 

n=—h{mod [d 2 ^l 2 ]) 
n=i;o(niod W) 




E 


"■ [dl.il] 

n=o'(mod [^ 2 ,^ 2 ]) 
n=/3'(mod W) 


2^{n) _|_ ^ Adi,d2Azi,Z2'''['^2, ^2 

{di ,lj ) — 1 


E 


2 ^{n) 


" [d 2 ,h] 

n=Q:"(mod [di,/i]) 
n=^" {mod W) 

(4.27) 


Here, we have taken out [di,li] and [^ 2 , ^ 2 ] respectively from the 1st and 2nd term in the last equality 
of 4.27 and replaced the summation over n + h with a summation over n. This is permitted because 
h is very small compared to x. Note that therefore a' is coprime to [^ 2 )^ 2 ]) cy" is coprime to 
and 13', (3" are both coprime to W. 


We assume that Arf^_d 2 ’s symmetric i.e \di,d 2 — ^d 2 ,di ^ 1 ,^ 2 . This will allow us to 

interchange the indices 1 and 2 in the second term in the last equality of (4.27) to get 


/ 


\ 


52 


E ^di 42 ^h,l 2 '''V^^dl] 

{di,lj)=l 


E 


2 ^{n) _|_ 


E 


(^{n) 


r , , nr^ r, 1 

n=a'(mod [^ 2 ,^ 2 ]) n=a"(mod [^ 2 ,^ 2 ]) 

\ n=^\mod W) n=^" {mod W) / 


We make use of Corollary 3.12 with m = VF[d 2 )^ 2 ] to evaluate the inner sums of the above 
expression. In order to apply the Corollary, the following conditions must hold 

/ \ 212,—e / 

Here, the second condition above holds because we have assumed to be symmetric. So, their 

support is also symmetric in the indices 1 and 2 . Since W is very small as compared to x, it can be 
swallowed into the term. So it is enough to have 

max{d^'^^(i 2 ; < -2 (4.28) 

where z = This explains why we have chosen the support of the sieve weights to be S(z). 



Now, Corollary 3.12 gives us 


C- _ 9 ‘F(^) sr^ \ \ '''[dldl] 1 

{di,lj)=l 


( 


logx-log[di,Zi]+c' + 2 E 


V 


p\W[d2,l2] 

p>2 


logP 

P-2 


+ Oe { X^/^+^ E 


\^di,d2 11 Azx,Z2 I 

lFV4[d,,/,]l/2[d2,Z2]V4 

di .Li 


(4.29) 
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where we have (with the notation of Corollary 3.12), 

p{p - 1) 


9{p) = „ 

p — 2 

pm 

where we have with /(n) = njT{n) that 


and c(W) = 


p^W 


{p - 1 )^ 

P{P - 2 ) 


We can therefore write 


S 2 = 2 x 


c{W){M 2 — M21 + M22) + E2 


,^2 42 


(4.30) 


^dl ,c /2 


M 21 = 




M 22 = 2 ^ 




log 


E 


logp 


(4.31) 


(di ,lj ) 


p>2 


E2 <€. y~i i'^'^i.c*2ii-^h42i 




Out of all these terms, only the terms M 2 and M 21 will be contributing to the main term. The rest 
of the terms M 22 and E 2 will be error terms. 

To compute these terms, we prove a few Lemmas concerning additive functions. The purpose of 
doing this is to simplify the computation for M 21 and the estimation for M 22 with these Lemmas. In 
both of these cases. Lemma 4.7 is directly applicable. In simpler words. Lemma 4.7 helps us give an 
asymptotic formula for M 21 and as well as we can use it to provide an upper bound for M 22 - 

Lemma 4.6. Let L{n) he an additive function defined on squarefree integers by L{n) = ^^L(p). Let 

p\n 

C'ri,r 2 as in (4.13). Then 


E 




f{di)g{d 2 )' 


di=0(mod Vi) 

where /i * 1 = / and gi* 1 = g. 
Proof. We have 

^di,d 2 


fi{ri)gi{r2) 


L(ri) + 2 ^ 


a 


pri,r 2 


Lip) 


fiiri)giir2) P - 2 


E 


di=0(mod ri) 


fidi)f{d 2 ) 


L(di) = ^ 


^di,d 2 


/(d,)9(<,,)Ei(p) = E^(i>) E 

dj=0(mod Vi) p\<il P 


^di,d 2 




Q 


lp,ri],r 2 


a 


ri,r 2 


di=0{mod [p,ri]) 
d 2 = 0 (mod r 2 ) 

c, 


fidi)gid 2 ) 


fi\p,ri]giir2) /i(ri)c/i(r2) 


L(ri) + 2 ^ 


pri,r 2 


Lip) 


h{ri)gi{r2) P - 2 


In the last step above, we have split the sum into two cases, depending on whether p | ri or not. This 
completes the proof. □ 


Proposition 4.7. Let Lin) be as in Lemma 4.6 with the further restriction that Lin) <C logn. Then 


^di,d 2 ^h,l 2 

^^^^^,f[difii]9[d2,l2] 


L[duli]=Yl 


^l ,'^2 


Ti 


fiiri)9iir2) 


L(ri) +4^ 


a 


ri,r 2 


Lip) ^ 


fiiri)9iir2) ^ P-2 

' l P 


pri,r 2 


+ 0 


log® 2 ; 
D 
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Proof. Now, 


LHS= 


,d 2 ,^2 


^^^^^J[duh]g[d2,l2] 


m,h]= Y1 


^di ,^2 ,^2 


M-1 
{_di,lj ) — J- 


■(L(di)+L(/i)-L(di,/i)) 


The fact that L[di,/i] = L{di) + L{li) — L{di,li) follows from the identity [di,li] = and the 

fact that L is additive. 


As we have done while computing Si, we want to get rid of the conditions (di,/2) = (^2,^1) = 1- So 
we multiply a factor 

Ksi 2 )g{s 2 l) 

Sl 2 \(<ii,l 2 ) 

S2l\{<i2,ll) 

We then have the same conditions as before that (risi2, ^2521), (?’iS2i, r2Si2) £ S{z). 

Moreover, one can write 

Y /i (^1)91 (’’2) 


f[di,li]g[d 2 ,l 2 ] fidi)f{h)g{d 2 )g{l 2 ) 


ri\{di,li) 


Therefore, 


Si " S fIJ \ n S fi(si2)fi(a2l) 

^h,h 


Y f‘(*12M»2l)/l('’l)9l(r2) Y 


E 


^di,d 2 


J{di)g{d 2 ) 

ti=0(mod riSji) j \ai=0(mod riSij) 


L{di) 


(4.32) 


There are two bracketed quantities in the last term above. The first one is Crj^s 2 i,r 2 si 2 divided by the 
quantity fi{riS 2 i)gi{r 2 Si 2 ). For the second bracketed quantity, we invoke Lemma 4.6. We therefore 
obtain the expression 


E. = E 


Ks 12 S 2 i) C, 


^ 1521,^2512 


{ 


/i(s12S2i)5i(s12S2i) /l (n)£/l (?’2) ^ 


L[riSi2)Crisi2,r2S2l + 2 E|^ 2^prisi2,r2S2i 

P ^ 


(4.33) 

Above, we have used that fact that fi{p) = (p —2)/2. Again, the contribution to (4.33) from all those 
(•512, S21) for which atleast one of Sij > D is 


« E 


1 


log^ 2 log® z 

< 


(si2,S2i)G5(2:) 

S12>D 

{ri,r 2 )eS{z) 


fi{si 2 S 2 i)giisi 2 S 2 i) fi{ri)gi{r 2 ) D 


where we have used the estimate Cri,r 2 log z from Proposition 4.4, the assumpion that L{n) <C log n 
and the fact that /i G II2, 51 G Hi. 

So, we can assume S12 = 1521 = 1 at the cost of error O ^ ^ ^. This gives 


us 


E, = E 


G 


ri,r2 




fi{ri)gi{r 2 ) 


Liri)Cr,,r2+2Y^Cprur2 ] +0 

p<z ^ 


D 
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By the same argument, one obtains 

We do similar operations as above for We have 

^di ,^2 '^^ 1 5^2 




piid.M 


Y1 E /i(si2)/x(s2l) 

^ij I {di ilj ) 

/ 


f^isi2S2i)fi{ri)gi{r2)'^L{p) 


Ti.Si^ 


(4.34) 


E 

dl=0(mod [p,risi2])) 
\ d2=0(mod r2S2i) 


A, 


d\,d2 


f{di)g{d2) 


E 


A 


h,l2 


1 Zl=0(mod [p,riS2l])) 
/ \ Z2=0(mod r2Si2) 


f{h)g{l2) 




Substituting the two bracketed quantities with appropriate C'ri,r- 2 ’s, we obtain 


E,= E 


h('Sl2'S2l) 


^3 r^.. fl{si2S2l)gi{si2S2l) ^ 




/l (^1)91 ^[p,riS 2 i],r 2 Si 2 

flip, nsi2]/i[p, riS2l]gi{r2S2l)gi{r2Sl2) 


(4.35) 


We now split the sum (4.35) into 4 cases depending on whether p | ri, p | S 12 , p\ S 21 or p f riSi 2 S 2 i- 
These are the only cases that can occur because we have p^(riSi 2 S 2 i) = 1. We would then have 


E,= E 


P(S12S21) 


C*risi2,r2S2iC*riS2i,r-2Si2 


"3 /l(si2S2l)Pl(si2S2l) /l(n)Pl(r'2) 

' ?, 1 7."? 


L(ri) 


_l_ 2 p('Si2S2i) L{p) C, 

riXj h{si 2 S 2 l)gi{si 2 S 2 l) ^ (p - 2) 


r\si2,r2S2l ^priS2i,r2S^2 


+ 2 E 


P(S12S21) 


E 


L(p) C, 


fi(ri)gi(r2) 

Cr 


prisi2,r2S2i'^riS2i,r2Si2 


fi{si2S2i)gi{si2S2i) ^ {p-2) h{ri)gi{r2) 


+ 4 E 


P(S12S21) 


/l(si2S2l)Pl(si2S2l) V (4’- 2)^ 


E 


L{p) 

^priSi2,r2S2l^priS2i,r2Si2 


fi{ri)gi{r2) 


(4.36) 


Of the four terms in (4.36) above, the 4th term is 0{\og^ z). This follows by using the estimates 
Cri,r 2 ^ log 2 ;, L{n) <C logn and the fact that fi G O 2 and gi G Oi. The 1st, 2nd and 3rd terms will 
be contributing to the main term. 

Again, the contribution to the 1st term, 2nd term and 3rd term of (4.36) when atleast one of 


Sij > D is 


< 


E 


1 


E 


log^ z 


< 


log® z 


, , fl{si 2 S 2 l)gi{si 2 S 2 l) . , /l(n)Pl(r- 2 ) " D 

{si 2 ,S 2 l)&S{z) {ri,r 2 )GS{z) 

S12>D 


So we may assume S 12 = '§21 = 1 at the cost of error O ^ ^ ^. Hence, 


E,, = E 


^2 

^ri,r2 


-L(ri) + 4^ 


C, 


ri,r2 


logp^ 


+ 0 


log® z 


'3 ^ fi{ri)gi{r2) /i(ri)pi(r2) ^p-2 \ D 

Putting together the expressions for E, E. and we get the desired result. 


□ 


We move forward with with our computations. We compute M 2 . 
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Evaluation of M2 

We multiply a factor of ^ Ai('Si 2 )A*(s 2 i) to get rid of the conditions (^ 1 ,^ 2 ) = {d 2 -,li) = 1- We 

I ih ) 

then have as before that /r^(rir 2 Si 2 S 2 i) = 1- Hence 


M 2 


E 


^di ,cZ2 


E 


A‘(s12'S2i) 


/i(s12'S2i)5i(s12S2i) 


E 

di 

E 




logx + c' /f-' f[di,h]g[d 2 ,l 2 ] f:^J{di)gid 2 ) f{h)g{l 2 ) ^ 

\dijtj) — 1 di,li '^i\\di^li) Sij\(^di,lj) 


n 


C'risi2,r2^riS2i,r2Si2 

fi{ri)gi{r 2 ) 


As before, the contribution to the above sum when atleast one of Sij > D is O ^ 
So we may assume S 12 = S 21 = 1 at the cost of error O ^ ^ ^. 


Therefore, 


M 2 


log X + c‘ 


7 = E 


r2 

ri,r 2 ^ Q 


fiiri)gi{r 2 ) 


log® X 
D 


E 


^ 2 (ri)/i 2 ('^ 2 ) / 2.27 x ,27 N/-, 2 /^logri log r 2 


(ri,r2)eS(z) 


/l(n) 51 (^’ 2 ) 


B^hi{n)hi{r2)Qi 


logz ’ log 2 


+ O(logz) ) + O 


(4.37) 


log® X 
D 


where we have substituted the expression for Cri,r 2 from Lemma 4.4. By applying Corollary 2.14 to 
the inner sum of the main term above, we obtain 


M2 = H®logx // siQl{si,S2)ds2dsi + 0 (log^ z)+ 0 


log® z 
D 


= (1 + o(l))B® logx JJ siQ|(si, S 2 ) ds2 dsi 
T 

where Q 2 is as given in (4.12). 


(4.38) 


Evaluation of M21 

Applying Proposition 4.7 to the expression for M 21 in (4.31) with L{n) = logn, we get 


fi2 

M 21 = ^ 


/i(n)5i(r2) 


log ri + 4 ^ 




ri,r2 


+0 


!\{r\)g\{r2) ^ p-2 


log®z 

D 


= M« + mS> + OC^ 


Substituting the expression for r 2 . we obtain 


1 > 7 ^ p‘^{r2) f ^2u2r \u2r ^ log n 3 log n logr 2 , , 

M2\'=logz 2^ TTrTTTTT ( ^ ^iAi)^2H2)-nWT<32 ( TWT’TWT ) +<^(log^) 

(ri,r2)GS(A -1^1 ^ 

By Corollary 2.14 applied to the above sum, we get 


log z ^ \ log z ’ log z 


= {1 + o{l))B^ logz JJ sIqI{si,S 2 ) ds 2 dsi 

T 


(4.39) 
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Next, we look at and substitute the value of C'pri,r 2 write Si = to get 


, , h{ri)gi{r2) , , , 

(r’l,r’2)G*(^) {pr^,r2)&S{z) 


^ E 


Ai(n)Ai(?’2)Cri, 


7’2 


(r’i,r2)G5(2) 


/i(n)5i(?’2) 


B E 


(pri,r2)&S{z) 


{p-l)‘^logp ( logp \ , r,n \ 

^-Q 2 Si + 1 -,S 2 + 0 (log 2 ;) 

p[p — 2) p \ log z 


(4.40) 


First, we focus on the inner sum of (4.40). We have the condition (pri,r 2 ) G S{z) which is equivalent 
to 

p < (p, hFrir2) = l 

(p - 1)^ 

with r]{s) as defined in (2.5). Observing that —--- = 1 + 0(l/p^), it follows that the inner sum 


of (4.40) is 


p{p - 2) 


« E 


‘“"-«.(*.+||f.«)+o(B x: i55i^i+o(iog.) 


p 


2<p<2’JG2)-'>i 

By Proposition 2.2, the main term above becomes 


p\Wrir2 


V{s2)-si 


V{S2) 


B{\ogz + 0{l)) j Q 2 isi+t,S 2 )dt = B{logz + 0{l)) j < 32 (^ 1 , 52 ) 


0 

ri{s2) i?(si) 


Si 


B(logz + 0(1)) j j P{ti,t 2 )dt 2 dti = B{logz + 0{l)) Qi{si,S 2 ) 


Plugging the expression for the inner sum back into (4.40), we have 

p2(ri)p2(^2) 


Mf=-4 j: 


Mri)gi{r2) 


{ri,r2)£S{z) 

X (B'^logz hi{ri)h2{r2) Qi 


log ri log r 2 
log z ’ log z 


/logn logr 2 \ 2 nA 

O G'li) 


+ o|iog=. V w 

The second error term in (4.41) turns out to be 0(log® z). This is done by replacing 

E by < ^ and ^ by ^ ^ ^ 

{ri,r2)GS{z) ri<z p\Wrir2 p\W p\ri p\r2 

r2<z 

and interchanging the order of summation. Coming back to the main term in (4.41), apply Corollary 
2.14 to the inner sum and get 


=-4{1 + o{l))B^ log z [ f siQi{si,S 2 )Q 2 {si,S 2 )ds 2 dsi 


( 4 . 42 ) 
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Estimation of M22 


For the estimation of M 22 , we can apply Proposition 4.7 with L{n) = - and obtain 

I p — 2 

p\n 

p >2 


M 22 = ^ 


c, 


ri,r2 


n 


fi{ri)gi{r2) 


/ 

/ 

\ 

C^ri,r2 

E 

logp 

p-2 

V 

p|ri 

\p>2 

/ 


\ 


+ 4 E 


logp 


2<p<z 


(p-2) 


2 Cpri,r2 


+ 0 


log®z 


Using the estimate Cri,r 2 logz in (4.43) we get 

( 


M.22 


log^ Z 


r\<z 

r2<z 


fi{ri)gi{r2) 




logp 

p-2 ' ^ ip - 2)2 

Pin ^ 2<p<^ ^ 

\p>2 


+ 0 


1^ 

D 


< log^ 2 ; I ^ 

log® z 


r2<z 


gi{r2) 


E 


logp 




P / 


+ log^ ^ + 


log® z 
D 


< 


D 


(4.43) 


(4.44) 


From (4.31), we have 


Estimation of E2 


E 2 <. \^dl,d2\\^ll,l2\ 




dili 


Using the estimate Xdi,d 2 ^ from Lemma 4.4 and writing [di, k] = ———, we get 


E 2 


a;l/2+e iQgS 2 


E 


(dj, li) 

(di, /i)^/^(d2, 


1UV4 , 1/2 ,1/4, 1 / 2 ,1/4 

(di,d2)e5(2) “1 “2 U ‘2 

{li,h)&Siz) 


Since <C (log log x)^/^ and log® z are small, these terms can be swallowed into the term x^. 

Therefore, 

(9i,g2)&s{z) 


{di,d2)€S{z) 

(h42)eS(2) 

(di 


( 51 , 92 ) 65 ( 2 ) 


E 


-l/2^-l/4 
'2 


Now, 




dxKzjgx 
\g2d2<zl{g\d\fl^ 


(4.45) 


/ 


- 1/2 


Y. d-Y^ Y E 

d\<zlg\ d\<zlg\ g2d2<z/{gidi)^/^ di<z/gi 

g2d2<zl{g\d\fl^ 

Y dr^ < 

d\<zlg\ 


z 1“/- 

52(51 
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Substituting the above expression into (4.45), we obtain 

gi<z 

92<zlg\^^ 


Under the choice ^ we get 

E 2 <C, (4 4g) 

Substituting 2 ; = and noting that c(W) = l + o(l), we get from (4.38), (4.39), (4.42), (4.44), 

(4.46) and (4.30) that 

O rp 

S 2 = -il + o{l))B^R 2 {P) (4.47) 

where 


R2{P)= / / Sl(3 + e- Si)Q 2 (si,S 2 )dS 2 dsi +4 / / SlQl(si,'S2)<32('Sl,'§ 2 ) ds2 c?si 


(4.48) 


T T 

Remark 4 . 8 . In the above expression (4.48) for R 2 {P), there is an e occuring. Since e can be made 
arbitrarily small, we do not consider it for our computations. 


This completes the Proof of Proposition 4.3. 


4.3 The Value of A 

Proposition 4.3 gives us the asymptotic expressions for Si and S 2 in terms of the symmetric differ¬ 
entiable (in each variable) function P. There are two more functions, namely Qi and Q 2 defined in 
terms of P (See 4.12). We make use of a computer program to optimise the choice of P to determine 
a suitable value of A. 

Let r](s) be as defined in (2.5). Recall that 


V{s) 


l-2s/3 ifs<3/5 

3/2(1-s) ifs>3/5 


(4.49) 


Moreover, note that for any integrable function T : T —>• M, one can write 


v{s2) riih) 


F{ti,t2) dt2 dti = // F{ti,t2) dt2 dti 


(4.50) 


Ti„ 


Using the above observations in (4.49) and (4.50) and using the expressions for Qi and Q 2 are in 
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terms of P, from (4.12), we obtain 


Qi(si,S2) = < 


I i-iii 


1-|S2 


P{ti,t2) dt2 dti + / / P{ti,t2)dt2dti 0 < Si, S2 < ^ 


Si S2 


3 S2 
5 


J j P{ti,t2)dt 


2dti 


SI S2 

1-|S2 |(l-tl) 


/ / P{h,t2) dt2 dti 


0 <Si<-<S 2 <l 

5 


0 <S 2 <-r<Si<l 

5 


SI S2 

L-|si 


Q2(s1,'S2) = < 


P{si,t2) dt2 0 < Si < 


S2 


ki-si) 


(4.52) 


P{si,t2)dt2 - < Si < 1 
5 


S2 


Applying the same observations (4.49), (4.50) to the expression for Ri{P) and R 2 {P) from Proposition 
4.3, we obtain the following (rather complicated) expressions 

2 


3 3/3 i_ 2 ,^ 

5 5 


Ri{P) = 


( f 


0 0 


1-|S2 f(l-ti) 





\Sl S2 


P{ti,t2) dt2 dti + I / P{tip2)dt2dti 


3 S2 

5 


ds2 dsi 


11-1*1 

/ f(l-si) l-|ti 

\ 

•// 

/ / P{ti,t2)dt2dti 

0 3 

5 

\ Si S2 

/ 


1 f(l-Sl) /l-fs2f(l-tl) 


+ 



ds2 dsi 


ds2 dsi 


(4.53) 


3 0 

5 


^ l-^Sl 


R 2 {P) = j j si{3-si) 


0 0 


3 3 
5 5 


0 0 




1 |(l-si) 




si(3 - si) 


P{si,t2)dt2 ds2dsi 


S2 


P{si,t2)dt2 P{ti,t2) dt2dti + I / P{ti,t2) dt2 dti j ds2dsi 


1 / 1 gSl 


+ 4 si 


0 3 

5 




1-|S2 f(l-tl) 


3 52 

5 




P{si,t2)dt2 j j P{tip2)dt2dti 


3 

1 5 


\ 


/l-|s 2 |(l-tl) 



/ P{sip2)dt2 


/ / P{ti,t2)dt2dti 

3 0 

5 

\ "2 

/ 


\ Si S2 

/ 


ds2 dsi 


ds2 dsi 


(4.54) 



















32 


Keeping in mind that P should be symmetric, we make the following choice for P. We set 

7 

P(.x, 2/) = ^ aij{x + y)\x^ + y^y ( 4 . 55 ) 

i,j=0 

The choice of the coefficients will optimised by means of a sage program. We evaluate Ri(P) and 
R 2 {P) in terms of these coefficients. We hrst convert the coefficients atj into variables of the form 
{ao, 01,02, ... ,063} via the bijective map (i,j) 8i + j and then integrate to obtain the expressions 
i?i and i? 2 - 


n=7 

N=(n+1)*2 
x=var("x") 
y=var("y") 
def pos(n,i,j): 

return (n+l)*i+j 

aa=list (var (’ a_7od ’ 7 i) for i in (0..(N)) ) 
p=0 

for i in range(0,n+l): 

for j in range(0,n+l): 

p=p+aa [pos (n, i, j)] * (x+y) "i* (x*2+y''2) * j 

tl=var(’tl’) 
t2=var(’t2’) 
sl=var(’si’) 
s2=var(’s2’) 
t=var(’t’) 


P(a,b)=p(x=a,y=b) 

Qll(si,s2)=integrate(integrate(P(tl,t2),t2,s2,l-2/3*tl),tl,si,3/5) 
+integrate(integrate(P(tl,t2),t2,s2,3/2-3/2*tl),tl,3/5,l-2/3*s2) 

Q12(si,s2)=integrate(integrate(P(tl,t2),t2,s2,l-2/3*tl),tl,si,3/2-3/2*s2) 

Q13(sl,s2)=integrate(integrate(P(tl,t2),t2,s2,3/2-3/2*tl),tl,si,l-2/3*s2) 

Rl=integrate(integrate(Qll(sl,s2) "'2,s2,0,3/5),si,0,3/5) 

+integrate(integrate(Q12(sl,s2)*2,s2,3/5,l-2/3*sl),si,0,3/5) 

+integrate(integrate(Q13(sl,s2)*2,s2,0,3/2-3/2*sl),si,3/5,1) 

Q21(si,s2)=integrate(P(sl,t2),t2,s2,l-2/3*sl) 

Q22(si,s2)=integrate(P(sl,t2),t2,s2,3/2-3/2*sl) 

Q41(sl,s2)=sl*(3-sl)*Q21(sl,s2)*2+4*sl*Qll(sl,s2)*Q21(sl,s2) 

Q42(sl,s2)=sl*(3-sl)*Q21(sl,s2)*2+4*sl*Q12(sl,s2)*Q21(sl,s2) 

Q43(sl,s2)=sl*(3-sl)*Q22(sl,s2)*2+4*sl*Q13(sl,s2)*q22(sl,s2) 


R2=integrate(integrate(Q41(si,s2),s2,0,3/5),si,0,3/5) 
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+iiitegrate(integrate(Q42(sl,s2),s2,3/5,l-2/3*sl),si,0,3/5) 

+integrate(integrate(Q43(sl,s2),s2,0,3/2-3/2*sl),si,3/5,1) 

R2=R2.expand 0 
R1=R1.expand 0 

In this code, Ri and R2 have expressions for Ri{P) and R2{P) stored in them. These expressions are 
actually quadratic forms in the coefficients oq, ai,..., oes- 


To show that (4.1) is positive, we must have A^i > 82 - By Proposition 4.3, we must have 


A > 


2R2{P) 


(4.56) 


R^{P) 

So, we need to minimize R2{P)/Ri{P)- This ratio is entirely dependent on the choice of function P 
and since it so happens with this choice of P (See 4.55) that both R2 and Ri are quadratic forms in 
the coefficients oq, oi,..., aga, we are essentially looking to minimize the ratio of two quadratic forms. 
This is a well known problem with the following solution. 


Theorem 4 . 9 . Let Ri = ae^Mia and R2 = a^M2a be two quadratic forms, where Mi and M2 are 
positive definite real symmetric matrices. Then the ratio R2/R1 is minimized when a is an eigenvector 
corresponding to the smallest eigenvalue of Mfi^M2. The value of the ratio at its minimum is this 
minimum eigenvalue. 


Proof. For a proof, see Lemma 7.3 [May, Pg 20]. 


□ 


We minimize this ratio using a sage program. We will let A to be the matrix corresponding to 
the form Ri and B corresponding to R 2 • We define them as follows 

, d^Ri , ^ d^R2 

= a a = a a 

oaidoj oaidoj 

Then it is clear that Aa = 2Ri and a^Ba = 2 R 2 . Moreover, the matrices A and B are clearly real 
symmetric matrices. If one looks at the expressions (4.2) and (4.3) for S'! and 82 respectively, it clear 
why both the matrices A and B (for Ri{P) and R 2 {P) respectively) are positive definite. 

So by Theorem 4.9, the the minimum value of R 2 {P)/Ri{P) is the smallest eigenvalue of C, which 
is computed below. 

A=matrix(QQ,N) 

B=matrix(QQ,N) 
for i in range(0,N): 

for j in range(0,N): 

A[i,j]=derivative(derivative(Rl,aa[i]),aa[j]) 

B[i,j]=derivative(derivative(R2,aa[i]),aa[j]) 

C=A.inverse()*B 
min(C.eigenvalues 0) 

6.290731135292344? 

Therefore, from (4.56), it follows that (4.1) is positive for any 

A > 12.5814622705847 (4.57) 


Remark 4 . 10 . It is of course possible to reduce the value of A somewhat further, if we take P{x,y) 
to be of a higher degree. But the improvement obtained will be very little and A is not likely to go 
below 12. Plus, the process can become increasingly time consuming for any computer program as 
we increase the degree of the polynomial P{x,y). 
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